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Abstract 

In (1), a result to algorithmically compute the topology types of the level curves of 
an algebraic surface, is given. From this result, here we derive applications based on 
level curves to determine some topological features of surfaces (reality, compactness, 
connectivity) and to the problem of plotting. 



1 Introduction 



The study of the level curves of an implicit algebraic surface (i.e. the sections 
of the surface with real planes parallel to the xy-plane) gives a clue on how 
the surface is like. Take for example the well-known example of the Whitney 
Umbrella, whose equation is x 2 — y 2 z = 0. It is clear that for z > the level 
curves consist of two intersecting lines; for z — 0, the level curve reduces to 
one line; and for z < 0, the level curves consist of just one real isolated point. 
From this simple analysis one may get a good mental picture of the surface. 
Moreover, it provides an initial analysis to the topology of the surface, which 
can be used later to provide topologically-correct plottings. For instance, in 
the case of the Whitney Umbrella one has that a topologically-correct plotting 
must show the so-called "handle" of the surface, i.e. the half-line corresponding 
to x = 0, y = 0, z < 0. Moreover, in addition to the shape of the level curves, 
one also has the ^-intervals corresponding to each topology type. 

The problem of computing the different topology types arising in the family 
of level curves of a given algebraic surface (together with the ^-values cor- 
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responding to each type) has been addressed, for instance in (|l|) and (jl2j). 
where algorithms to compute these topology types are given. Moreover, in (|l|) 
applications of this approach to the study of offset curves are also given; see 
also (0) for further research on this topic. Now, in this paper, we apply these 
results on level curves to determine topological information on the surface, 
which may be of help for computing reliable plottings from the topological 
point of view. 

In this sense, the first application that we consider is to algorithmically de- 
cide whether a given implicit algebraic surface is real. That is, we provide 
an algorithm to check whether the intersection of the surface with M 3 is a 
two-dimensional set (in the Euclidean topology) or not. Furthermore, in the 
negative case, the algorithm can also be used to analyze whether the real part 
of the surface is empty (like for example x 2 + y 2 + z 2 + 1 = 0), or it consists 
of finitely many points (the case of x 2 + y 2 + z 2 = 0, whose real part is the 
origin) or it is a space curve (the case of x 2 + y 2 = 0). This kind of surfaces 
whose real part is not 2-dimensional may get error messages when one tries 
to draw a plotting. 

The second application concerns to the compactness of the surface. In this 
case, since one works with implicit algebraic surfaces, which are closed over 
the usual Euclidean topology, the algorithm essentially checks whether the 
surface is bounded w.r.t. the variables x, y, z, respectively. 

The third application directly concerns to the problem of plotting. In order 
to draw a plotting of a surface, the user has to introduce as an input a "box" 
/ = [ai, a 2 ] x [bi, 62] x [ci, C2] C IR 3 , so the output shows the part of the surface 
lying inside the box. Now, if the user is interested in computing a plotting 
where the main topological features of the surface are shown (i.e. which makes 
clear how the surface is like), some previous information must be known in 
order to properly choose /. Using the information on the topology of the level 
curves of the surface w.r.t. the variables x,y,z, we provide an algorithm to 
compute a box of this kind. 

The fourth application is a symbolic-numeric algorithm to compute the con- 
nected components of a surface by using level curves. Here, the algorithm 
determines how the level curves are connected to each other, and as a con- 
sequence also provides the number of connected components of the surface. 
Since the algorithm in fact computes how to join sections of the surface, we 
think that it might be used for plotting purposes. This question might be 
explored in a future work. 

The structure of the paper is the following. The first section contains some 
preliminary notions and related results on level curves. The second section is 
devoted to the problem of checking whether a given implicit algebraic surface 
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is real, or not. The third section analyzes compactness. The problem of com- 
puting a suitable box for plotting a surface (so the output shows the more 
relevant topological features of the surface) is addressed in the fourth section. 
The fifth section concerns to connectivity. 



2 Preliminaries on Level Curves 



In this section we provide some preliminary notions and results concerning 
to level curves, that are taken from (jl|). Thus, we refer the interested reader 



to (jl|) for further information. Moreover, here we also fix the notation to be 
used along the paper, together with some hypotheses to be requested on the 
surface S. 

In the sequel, we consider an algebraic surface S defined by an square-free 
polynomial F e M[x, y, z] with no univariate factor only depending on the 
variable z, i.e. 5" has no component which is a plane parallel to the xy-plane. 
In this situation, the level curves of S are the (plane) curves which are obtained 
by intersecting S with planes normal to the z-axis. Furthermore, given ieK 
we will denote the level curve corresponding to the plane z — b = by Sb- 

Note that one might similarly define the level curves corresponding to the 
x-axis and the y-eads, respectively. Thus, when necessary we will speak of 
£-level curves, where £ G {x, y, z}. Note that the topology type of a £-level 
curve can be described by computing a graph homeomorphic to it, which is a 
well-studied problem (see for example (jl), (0) and many others). Essentially, 
the vertices of this graph are the real critical points of the curve (i.e. the real 
singular points and the real points where the tangent is vertical), and the 
edges correspond to the branches of the curve joining two vertices. 



From Hardt's Semi- Algebraic Triviality Theorem (see (j^)) it can be derived 
that the number of topology types of the level curves of S, is always finite. 
In case that S is compact and non-singular the problem of determining these 
topology types can be solved by using Morse Theory (see(jl), (fill)). In the 
more general case of singular surfaces, two approaches can be considered. 
The first one comes from Differential Topology and uses elements of Whitney 



Stratification Theory (see p)). This approach has been used in (11 21 ). The 
second one comes from Computer Algebra and uses as an essential tool the 
notion of delineability (see (0)). This second approach has been developed in 
([]]). In the rest of the section, we recall some notions and results in (jl|) related 
to the computation, by means of (jl|), of the topology types of the level curves 
of S. 

Definition 1 We say that a el is a Critical Level Value if the topology of the 
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level curves of S changes at z — a, i.e. V e > there exists a e G (a — e, a + e) 
s?/c/j t/iat #ie level curves corresponding to a and a t , have different topology 
types. Moreover, we say that A C R is a Critical Level Set of S, if A is finite 
and it contains all the critical level values of S. 

Note that since the number of topology types of the level curves of S is always 
finite, the number of critical level values is also finite and therefore a critical 
level set always exists. Moreover, once a critical level set A has been computed, 
the topology types of the level curves can be obtained. Indeed, writing A = 
{o!i, . . . , a r }, we can decompose the z-axis as 

(—00, ai) U {«i} U a 2 ) U ■ • ■ U (a r _x, at r ) U {a r } U (a r , 00). 

Thus, taking a z-value for each open interval of the partition, and applying the 
existing algorithms for computing the topology type of a plane algebraic curve 
(see (ji|), (0)), one determines the topology type for all the z- values correspond- 
ing to the interval. The remaining finitely many level curves, corresponding 
to F(x, y, cxi), i = 1, . . . , r, are also analyzed with the same strategy. 

Therefore, the problem of determining the topology types of the level curves 
of S reduces to the computation of a critical level set. Let us recall the results 
in (jl|) which allow to do this. For this purpose, throughout this paper, besides 
of the hypotheses required above (i.e. S is implicitly defined by an square-free 
polynomial F G M[x, y, z] having no univariate factor only depending on z), 
we also assume that the leading coefficient of F w.r.t. y does not depend on 
the variable x. Note that this requirement can always be fulfilled by applying 
if necessary a rotation around the z-axis, which does not modify the topology 
of the level curves of S. 

Now, we consider the following notation: D W (G) denotes the discriminant of 
a polynomial G w.r.t. the variable w, i.e. D W {G) = Res^(G, |^), y/G denotes 
the square-free part of a polynomial G, and: 



Then, the following result holds (see (Jl|)): 

Theorem 2 Let S satisfy the above hypotheses. Then, it holds that: 

(1) If R(z) is not identically zero, then the set of real roots of R{z) is a critical 
level set. 





otherwise 








4 



(2) If R(z) and M are identically zero, then the set of real roots of D X (F) is a 
critical level set. 

(3) If R(z) is identically zero but M is not identically 0, then the set of real 
roots of M{z) is a critical level set. 

Remark 1 // R(z) is a non-zero constant, then there is just one topology 
type for all the level curves. Similarly for the case when M(z) is a non-zero 
constant. 



3 First Application: Reality of Algebraic Surfaces 



In this section we deal with the problem of deciding whether an algebraic 
surface is real of not; i.e. whether S H M 3 has dimension two with the usual 
Euclidean topology. For this purpose, we will show that the problem of check- 
ing the reality of S can be reduced to checking whether finitely many level 
curves of 5* are real, so the problem in R 3 can be reduced to a problem in 
M. 2 . Observe that in order to check whether an alg ebraic curve is real, one 



can use well-known algorithms, like for example (1141 ) . The following theorem, 
which can be found in (0), is essential. Here, we consider the usual definition 
of regular point of a surface F(x,y,z) = 0, i.e. a point P of the surface is 
regular if some first partial derivative of F P does not vanish at P. 

Theorem 3 S is a real surface if and only if it has at least one real regular 
point. 



Proof. See Theorem XI. 3. 6 in (J8|). 

From this theorem, one may derive another result, concerning to the level 
curves of S, which can be used to algorithmically check whether S is real. In 
order to see this, we need the following previous result. We denote the partial 
derivatives of F w.r.t. the variables x,y,z, respectively, by F x , F y , F z . 

Proposition 4 Let F e M[x, y, z], with deg y (F) > 0, be the defining poly- 
nomial of the surface S , and let A be the critical level set of S computed by 
means of Theorem If a Gl verifies some of the following conditions: 

(i) a is a root of the leading coefficient of F w.r.t. y, 

(ii) a is a root of the leading coefficient of M(x,z) w.r.t. x, 
(Hi) the polynomial F (x , y , a) has multiple factors, 

then a G A. 

Proof. In order to prove the result, we distinguish the cases when the poly- 
nomial R, defined in Section [21 is identically 0, or not. Let us see first the 
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case when R ^ 0. Here, we observe that if a is a root of the leading coefficient 
of F w.r.t. y, then by using the Sylvester form of the resultant Kes y (F, F y ), 
one has that z — a is a factor of D y (F) = Kes y (F,F y ) and therefore of M. 
Thus, it is also a factor of the leading coefficient of M w.r.t. x, and conse- 
quently, (i) implies (ii). Now, let us see that the condition (ii) implies that 
R(a) = 0, and therefore a G A. Indeed, using again the Sylvester form of the 
resultant Res x (M, M x ) = D X (M) = R(z), one deduces that if z — a is a factor 
of the leading coefficient of M w.r.t. x, then it is also a factor of R(z). Thus, 
R(a) = and a G A. Finally, let us see that if a verifies (iii) but not (i), then 



R(a) = 0. Let H(x,z) = D y (F); then M(x,z) = JH(x,z). Since a does not 
verify (i), the discriminant of F w.r.t. y behaves properly under specializations 
when z = a, i.e. D y (F a ) = H(x,a) (see (1151) ). However, since F a has multiple 
factors, we get that D y (F a ) = 0, so H (x, a) = 0. Therefore, z — a is a factor of 
D y (F), and consequently (ii) occurs. Hence, in this case, the result is proved. 
Now, let us see the case when R = 0. Since deg^(-F) > 0, then by Theorem 
Hit holds that M/fl but M G R[z]. Hence, by Theorem El A is the set of 
real roots of M, so clearly if a satisfies (ii), then a G A. Thus, let us see (i). 
For this purpose, let A(z) be the leading coefficient of F w.r.t. y, and let a 
be a real root of A. If deg y (F) > 1, then reasoning like before we get that 
z — a is a factor of M, and hence a G A. If deg y (F) = 1, then |^ = A(z), 
and ReSj / (F, ^) = A(z). Hence, in this case it also holds that z — a is a factor 
of ReSj / (F, and consequently of M. Therefore, we also get that a G A In 
order to prove (iii), we proceed as in the case R ^ 0. □ 



Now, we can prove the following result concerning to the level curves of S. 



Theorem 5 Let A be a critical set of the surface S determined by applying 
Theorem d Then, S is real if and only if there exists at least one real level 
curve S a of S, with a G M and a ^ A. 

Proof: If S is real, then, by Theorem [31 there exists a regular real point 
P G S. Thus, the implication (=^>) follows from Implicit Function Theorem. 
Let us consider the implication (<=)■ For this purpose, we separately analyze 
the cases when F G ~R[x, z], and when F depends on the variable y. We start 
with the case F G M[x, z\. In this situation, let C xz be the plane algebraic 
curve defined by F in the xz-plane. Thus, S is real iff C xz is real. So, in order 
to prove that S is real, it suffices to prove that C xz is. Let us see that this 
holds. Now, by hypothesis there exists a G R such that a £ A, and verifying 
that the corresponding level curve S a is real. Since F does not depend on 
y, the level curves of S are lines normal to the xz-plane. Therefore, since S a 
is a real curve, one has that the intersection point of S a with the xz-plane, 
which we denote also real. By Theorem [21 A is the set of real roots 

of the discriminant D X (F). Thus, a is not a root of D X (F). Therefore, P a is 
not a singular point of C xz and consequently C xz is a real curve. Thus, {<=) 
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holds for the case when F = F(x,z). Finally, let us see that (<=) also holds 
when F depends on the variable y, i.e. when deg^(-F) > 0. In order to see 
this, let S a be a level curve of S, real, and corresponding to the intersection 
of S with the real plane z = a, where a ^ A. Since a £ A, by Proposition 
H] the polynomial F a (x,y) = F(x,y,a) is square-free. Thus, since S a is real 
and F a (x, y) is square-free, we have that S a has at least one real non-singular 
point (x a , y a ) £ S a . Then, (x a , y a , a) is a real non-singular point of the surface 
S (note that V(F a ) = (F x (x,y,a),F y (x,y,a))), so by Theorem [3] the surface 
S is real. □ 

This theorem can be used to derive an algorithm for checking the reality of 
an algebraic surface. For this purpose, note that the condition in Theorem [5j 
i.e. the existence of a real level curve of S corresponding to a non-critical level 
value, can be tested by checking the reality of the level curves corresponding 
to intermediate z-values in between two consecutive critical level values. More 
precisely, one has the following algorithm: 

Algorithm: (Reality of an algebraic surface S) Given an algebraic surface S im- 
plicitly defined by a real polynomial F(x, y, z), square-free, with no factor only 
depending on the variable z, and such that lcoeff 2/ (F) does not depend on the 
variable x, the algorithm decides whether S is real. 

(1) Compute a critical set of S, A = {ai, . . . , a r } by means of Theorem [3 Let 
clq = — co, a r+ i = co. 

(2) Check whether there exists % £ {0,...,r} such that the plane algebraic 
curve defined by F(x,y,£i), where £j is taken in the interval (a^dj+i), is 
real. If it is, then return <C S is real^> else return <C S is not reals>. 

Remark 2 Note also that in case that the surface is not real there would be 
three alternatives: (i) the real part of the surface reduces to a space curve; 
(ii) it consists of finitely many points; (Hi) it is empty. Moreover, one can 
algorithmically decide which is the case by inspecting the level curves. More 
precisely, in case (Hi) all the level curves are empty; in case (ii), there are 
just finitely many non-empty level curves, all of them corresponding to z- 
critical level values, and consisting of finitely many real points. Finally, case 
(i) is identified when (ii) and (Hi) do not happen, and all the level curves 
corresponding to non-critical z-values are either empty or consisting in finitely 
many real points. 

Example 1 Let S be the algebraic surface defined by 

F(x, y, z) = (x 2 - l) 2 + (y 2 - l) 2 + (z 2 - l) 2 - 3/2. 

Note that S satisfies the imposed hypotheses. Let us see whether S is real. For 
this purpose, we apply Theorem^ to obtain the following z-critical level set of 
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S: 

A z = {-1.491557867, -1.306562965, -0.5411961001, 
0.5411961001, 1.306562965, 1.491557867} 

Now, we check if there exists some level curve, corresponding to a z-value not 
in A z , which is real. For z < —1.491557867 we get that the level curves are 
empty overM. 2 , but for z = —7/5, which is intermediate between —1.491557867 
and —1.306562965, we get the z-slice 

{(x 2 - l) 2 + (y 2 - l) 2 - 723/1250 = 0, z = -7/5} 
which is real. Therefore, we conclude that S is real (see Figure 1). 
Example 2 Consider the surface S defined by 

F(x, y, z) = x 4 + 2x 2 y 2 - 2x 2 + y* - 2y 2 + 1 + z 2 . 

Note that S satisfies the imposed hypotheses. In this case, a z-critical set of S 
is A z = {0}, i.e. the z-slices of S have at most three different topology types. 
However, the z-slices for z = — 1 and z = 1 are empty curves over IR, so from 
Theorem [H we deduce that for z > and z < the surface is empty over the 
reals. Therefore, S is not real. In fact, the only real points of the surface are 
the points of the z-slice corresponding to z = 0, which is the circle 

{(x 2 + y 2 -l) 2 = 0,z = 0} 



4 Second application: Compactness 

Here we show how to use level curves to algorithmically decide whether S is 
compact. Now, since S is implicitly defined by a polynomial F e M.[x, y, z], 
then it is obviously closed. Thus, in order to check whether it is compact, it 
suffices to check whether it is bounded, which is equivalent to decide whether 
it is bounded w.r.t. the x, y and z variables, respectively. 

Let us see how to check whether S is bounded w.r.t. the variable z. For this 
purpose, let A = {ai, . . . , a r } be a critical level set of S, where a\ < ■ ■ ■ < a r . 
Then, 5* is bounded w.r.t. the variable z iff for z > a r and z < ai, the level 
curves of 5* are empty over M. 2 . Moreover, since by Theorem [2] the topology 
type of the level curves of 5" stays invariant for z > a r , and also for z < a\, 
in order to check whether the condition holds it suffices to take Zq < a\ and 
z r+ i > a r , and then to analyze whether the level curves S Zo , S Zr+1 are empty 
or not over M 2 . Note that for this purpose one may adapt the strategy for 
deciding whether a given algebraic curve is real. 
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Similarly for the x and y variables. However, observe that in order to compute 
^-critical sets, with £ e {x,y,z}, by means of Theorem [21 one needs that the 
hypotheses of Theorem [2] hold not only for the variable z, but also for x, y, 
respectively. To ensure that this happens, one may always apply if necessary a 
linear transformation so that the polynomial F e M[x, y, z] defining S has no 
univariate factor, and verifies that \coeff x (F), lcoeff y (F) and lcoeff z (-F) are all 
constant. Observe that this kind of transformations preserves the topological 
properties of the surface. 

Thus, one may derive the following algorithm: 

Algorithm: (Compactness of an algebraic surface S) Given an algebraic surface 
S implicitly defined by a real polynomial F(x,y, z), square- free, with no uni- 
variate factor, and such that lcoeS x (F) ,lcoeS y (F) ,\coeS Z (F) are constant, the 
algorithm decides whether S is compact. 

(1) Compute an upper bound k z of the absolute value of the elements of a 
critical set of S. 

(2) Check if the plane algebraic curves defined by F(x, y, —k z — l) and F(x, y, k z + 
1) are both empty over M 2 . If this does not happen, then return <C S is 
not compacts. 

(3) Proceed in an analogous way for the variables x and y. If all the tested plane 
curves are empty over M 2 , return <C S is compacts, else return <C S is 
not compacts. 

Remark 3 In this case one does not need to compute the real roots of the 
polynomial provided by Theorem [H, but just upper and lower bounds on them. 
This can be done by applying existing algorithms (see for example 

Example 3 Let S be the surface in Example [21 which fulfills all the require- 
ments of the algorithm before, and let us see whether it is compact. For this 
purpose, we have to check whether it is bounded. This is equivalent to checking 
whether the z-slices below the least z-critical level value and above the greatest 
z-critical level value are both empty over M 2 ; similarly for y and x. In this 
sense, for z = —2 and z = 2 one gets 

{(x 2 - l) 2 + (y 2 - l) 2 + 15/2 = 0,z = -2} 
{(x 2 - l) 2 + (y 2 - l) 2 + 15/2 = 0, z = 2} 



which are obviously empty overM 2 , 
it also holds for y and x. Therefore, 
it is compact (see Figure 1). 



so the condition holds for z. By symmetry 
we deduce that S is bounded, and therefore 
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5 Third application: Plotting Boxes 

Here we address the problem of computing an interval / = [—a, a] x [—b, b] x 
[— c, c] cK 3 , so that the plotting of S in / shows the main relevant topological 
features of 5*. For this purpose, the information on the £-level curves of S, 
£ G {x, y, z}, is used. 

More precisely, we consider the following definition, which provides a criterion 
to compute /. 

Definition 6 We say that the interval [—m x ,m x ] x [—m y ,m y ] x [—m z ,m z ] C 
M 3 is suitable for plotting S if, for £ G {x,y,z}, —m^,m^ are not ^-critical 
level values of S and [—m^,m^\ contains all the ^-critical level values of S . 

Thus, if I is "suitable for plotting" S, one can be sure that out of / there 
is no change in the topology type of the £-level curves of S\ Note that the 
computation of a suitable I requires to compute critical level sets for the 
variables x, y, z, respectively, so one requests the same hypotheses as in Section 
HJ Observe also that Remark [3] also holds for this case. Thus, the following 
algorithm is derived: 

Algorithm: (Suitable interval for plotting an algebraic surface S) Given an alge- 
braic surface S implicitly defined by a real polynomial F(x,y,z), square-free, 
with no univariate factor, and such that lcoeff^-F), \coe& y (F), lcoeff 2 (.F) are 
constant, the algorithm determines a suitable interval /cl 3 for plotting S. 

(1) For £ G {x, y, z} compute an upper bound of the absolute values of the 
elements of a ^-critical set. 

(2) Return the interval / = [— k x — 1, k x + l] x [— k y — 1, k y + l]x[—k z — l,k z + l]. 

Example 4 Consider again the surface in Example^ Here, we have that 

A z = {-1.491557867,-1.306562965,-0.5411961001, 
0.5411961001, 1.306562965, 1.491557867} 

is a z-critical level set of S. Furthermore, by symmetry, 

A z A x 

Thus, the interval 

I = [-1.5, 1.5] x [-1.5, 1.5] x [-1.5, 1.5] 

is suitable for plotting S. The picture of the part of S lying in I is shown in 
Figure 1. 
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Fig. 1. F(x, y, z) = (x 2 - l) 2 + (y 2 - I) 2 + (z 2 - I) 2 - 3/2 
6 Fourth application: Connectivity 

There exist purely symbolic algorithms, based for instance on the notion of 
road map, to compute the connected components of S (see (jH) for further de- 
tails. Here we propose an alternative approach to solve this question by means 
of a symbolic-numeric algorithm based on level curves. The empirical perfor- 
mance of this new method is quite satisfactory. The main idea of the algorithm 
is to determine how the connected components of a level curve corresponding 
to a non-critical z-value join to the connected components of the level curves 
corresponding to the critical level values immediately below and above, re- 
spectively. For this purpose, we take points on each connected component of 
the non-critical z-value, and then we generate a space curve (as the solution 
of a system of differential equations) which connects it with some connected 
component of the critical z-level curve immediately below/above (see Figure 
2). Thus, once we know how to join the connected components of the level 
curves corresponding to non-critical and critical ^-values, the number of con- 
nected components of S can be obtained as the number of "connected chains" 
(whose elements are connected components of level curves) computed in the 
process. 

In order to solve this problem, we require some more conditions on the surface 
S to be analyzed. More precisely, the following hypotheses must be satisfied: 

(i) S is defined by an square-free polynomial F, with no factor just depending 
on the variables y, z. If this holds, then gcd(F, F x ) = 1, so the variety C 
defined by F = F x = is a space algebraic curve; recall that F x denotes the 
partial derivative of F w.r.t. x. 

(ii) There does not exist any plane z — a = containing infinitely many points 
of the curve C introduced in (i), i.e. not containing infinitely many points 



of S where F x vanishes, 
(iii) S is not asymptotic to any plane of equation z — b = 0, where b is a critical 
level value. 

Observe that in case that S does not fulfill some of these conditions, almost 
all linear transformations lead to a new surface where the three requirements 
hold. Note that linear transformations preserve the topological features of the 
surface. 

Now, in the sequel let A = {ai, . . . , a r }, where a± < ■ ■ ■ < a r , be a z-critical 
level set of S; furthermore, we set a = — oo, a r+ ± = +00. Moreover, let 
b < ■ ■ ■ < b r verify a« < bi < ai + \ for all i G {0, . . . , r}. With this notation, 
our problem is to decide, for each ai,bi,a i+ i, how to connect the connected 
components of the level curves S^.S^, and Sb v S ai+1 , respectively. Observe 
that computing the topology graph of a level curve (which is planar) one 
obtains the connected components of it (see for instance (0)). Also, note that 
it may happen that the level curves for values in some of the intervals (a,, bi) 
are empty over the reals, in which case there is no need of connecting S^. , S ai , 
and , S ai+1 . In addition, observe that, since S verifies condition (iii), every 
connected component of some level curve of S corresponding to a non-critical 
z-value, joins to some connected component of the level curve corresponding 
to the ^-critical level value immediately below (resp. above). In fact, this is 
the reason why we request condition (iii). 

For this purpose, the strategy is to use a symbolic-numeric algorithm which 
essentially works as follows: 

Algorithm: (Connected components of an algebraic surface S) 
Given an algebraic surface S implicitly defined by a real polynomial F(x,y,z) 
verifying the conditions (i), (ii), (iii), the algorithm determines the number of 
connected components of S and a description of them in terms of level curves. 

(0) Compute the topology graph of (if S& 4 fllR 3 = take another i), and the 
singular points of S ai ; the information on the singular points of S ai will be 
used at step (2), in some cases (see Remark H]), and also at step (3)). 

(1) Take a real point P in each connected component of Sb r 

(2) Use a path continuation method to connect P with some point Q in S ai , to 
be computed by the algorithm; in order to do this, we travel from P to Q 
by following a space curve, contained in the surface, which is computed as 
the solution of a system of differential equations. 

(3) Identify the connected component of S ai where the final point Q belongs 
to. For this purpose, we compute the topology graph of S ai introducing the 
point Q as a vertex of the graph. 

(4) Join by an edge the starting connected component of and the finally 
reached connected component of S ai . 
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(5) Proceed in an analogous way to connect P with some point Q* in S ai+1 . 

(6) After carrying out the computation for all the connected components of all 
the Sfy, one gets several connected chains, each one consisting of some con- 
nected components of the S^, S ai joined by edges. The number of connected 
chains, is the number of connected components of the surface. 

Observe that in the execution of step (0), points on each connected component 
of Sbi are computed (see e.g. (0) or (0)), so step (1) can be executed afterwards. 
Now, let us describe with more detail step (2). We consider the solution to 
the following system of differential equations: 

( F 

fx 

f y' = F x 

z' = -1 
^ x(0) = Xi) y(0) = yf, z(0) = bi 

where P := (xj, y^ bi) G S^. In case that this differential system has a symbolic 
solution (x(t),y(t), z(t)), one may see that it corresponds to a space curve 
contained in S. Moreover, since z'(t) = — 1 and 5* fulfills condition (iii), then 
this space curve reaches z = a^. Furthermore, since one may decompose the 
part of S with (eij, bi) into non-intersecting pieces, each one corresponding 
to a different connected component of (see (jH) for a careful proof of this 
fact), the choice of the initial point for a particular connected component of 
does not affect the connected component finally reached. In other words, 
the connected component reached at z = ai is always the same for all the 
points of a same connected component of S^ (see also Figure 2). 



Starting point 




Final point 



Fig. 2. Idea of the connectivity algorithm 
However, in general the differential system above may not have a symbolic 
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solution, so numerical methods must be applied. In our case, we used the 
package of maple for numerically integrating differential equations. Here, one 
may see that as the numerical method goes on, the error considerably grows, 
so the point of z = finally reached cannot be recognized as belonging to 
any connected component of S ai . For this purpose, at each step the solution 
provided by numerical integration must be corrected. In order to do this, each 
solution P^ k = (x i>k , y^k, z i>k ) is corrected to P^ k = (x iik , y~i,k, Zi,k) by computing 
a point of the level curve h close to the solution (xi ;k , y^k, z^h) (observe that 
the ^-coordinate is the same in P i k and in P^k)- For this purpose, we take the 
line passing through P ik in the direction of v = (F x (P i k ), F y (P ijk )), and we 
compute the intersection points of this line with k ■ This new point is used to 
go on with the numerical integration process. Although we have not analyzed 
the convergence of the method, the empirical results we get in this way are 
satisfactory. In this sense, the following remark must be taken into account. 

Remark 4 If K is a singular point of S ai ( notice that these points are com- 
puted in the initial step of the algorithm), then we assume that it is reached 
whenever the distance between P ik and K is smaller than a sufficiently small e 
previously fixed. Thus, in this case the computation stops and we assume that 
the starting point of is connected with K, i.e. that Q = K. 

In addition, there are two more situations which must be examined carefully: 

• If, before reaching the level plane z = a^, the numerical integration process 
hits a point where F x vanishes, then the method fails. This situation can be 
prevented by detecting whether \F x (x^ k , y^, Zi t k)\ < e; if this happens, we 
choose a different point of Ss i k to go on with the numerical process. 

• It may happen that S contains a 1-dimensional subset C of singular points, 
where £ is "isolated" in the following sense: given any point P G C, there 
exists a Euclidean neighborhood E p of P such that every point of S'nE'pniR 3 
is also a point of C For example, the handle of the Whitney Umbrella 
x 2 — y 2 z = 0, which is obtained for negative values of z, provides an example 
of this situation. Unless C is parallel to the xy-plane, this phenomenon can 
be detected by identifying the presence of isolated points in non-critical level 
curves. 

Hence, assume that has some isolated point B, therefore belonging 
to some C The points of C are singular points of S, so F x vanishes at B 
and therefore the system of differential equations before cannot be used 
to compute the point A in S ai which must be connected with B. Now, in 
this case we use the projections of the curve C onto two coordinates planes 
to compute how the points of C in z = bi and z = cu, respectively, are 
connected, in analogy with the method described in (|2|) to compute the 
topology of space algebraic curves. More precisely, we consider a coordinate 
plane so C is not normal to it. We project C onto this plane, we determine 
the projection tt(B) of the point B, and then we determine the projection 
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tt{A) of the point A by either applying a path continuation method (over 
the projection), or by computing the part of the graph associated with the 
projection between z = hi and z — dj. Carrying out this process on two 
coordinate planes, the point A is obtained. 

Finally, in order to connect z = hi and z = czj+i, we apply an analogous process 
to the differential equation system: 

( F 

X= - Fy -F~ x 

f y' = F x 

z' = 1 

^(0) =2*; 2/(0) =yi',z(0) = hi 

Note here that the third equation is different from the system before (z' = 1 
instead of z' = —1), since in this case one has to move "up" from z = hi to 
z = a i+ i. One may see that also the first equation has changed. However, the 
space curve that one obtains by integrating these equations lies also in the 
surface S. 

Example 5 Consider the algebraic surface S defined by F(x, y, z) = x 2 + 
y 2 + z 2 + 2xyz — 1. A z-critical level set of S is A z = {—1, 1}. Because of 
the symmetry of the surface, we have that A x = A y = A z , so for example 
[—2, 2] x [—2, 2] x [—2, 2] is suitable for plotting S. A plotting of S in this 
interval can be seen in Figure 3; this figure was computed with maple. 



■ ■ 

Fig. 3. The cubic x 2 + y 2 + z 2 + 2xyz — 1 = 
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Observe that, from Figure 3, it is not completely clear whether the surface is 
connected or not. However, using the algorithm we check that S is connected. 
Indeed, Figure 4 shows the different topology types corresponding to the z-level 
curves for the cases z < — 1, z = — 1, —1<z<1,z = 1,z>1, respectively. 
Moreover, also in Figure 4, each connected component of each z-level curve has 
been joined with the corresponding connected component of the z-critical level 
curve immediately above/below, according to the algorithm in this section. The 
connections between connected components of level curves are represented by 
dotted lines. Here, one may see that there is just one connected "chain", formed 
by the connected components of the level curves which are joined one another. 
Hence, S has just one connected component, and therefore it is connected. 





Fig. 4. Connections between the level curves of x 2 + y 2 + z 2 + 2xyz — 1 = 

Example 6 Consider the quartic (x 2 — l) 2 + (y 2 — l) 2 + (z 2 — l) 2 — 3/4 = 0. 

In this case, a z-critical level set of S is 

A z = {-1/2 - y/3/2, 1/2 - y/3/2, -1/2 + y/3/2, 1/2 + y/3/2} 



The different topology types for the z-level curves are shown in Figure 5; more- 
over, one may check that for z > 1/2 + y/3/2 and z < —1/2 — y/3/2, the level 
curves are empty over M. 2 . Furthermore, also in this picture we have repre- 
sented (in dotted lines) how the connected components of these level curves 
are joined with each other. This information has been computed by using the 
algorithm in this section. 
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Fig. 5. Connections between the level curves of (x 2 — l) 2 +(y 2 — l) 2 +(z 2 — l) 2 — 3/4 = 

Since in Figure 5 one may see 8 different connected chains of components of 
level curves, one has that S has 8 connected components. A plotting of this 
quartic can be seen in Figure 6. 




Fig. 6. The quartic (x 2 - l) 2 + (y 2 - l) 2 + (z 2 - l) 2 -3/4 = 
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